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The internal friction in the crystals grown in two modes, namely, slow and anomalously fast 
(burst-like growth) modes, is measured in the temperature range 0.49-0.75 K at a frequency of 
about 75 kHz. An additional contribution to the damping decrement and softening of the dynamic 
modulus are detected, and their relaxation to equilibrium values at a time constant of about 3 ms 
is observed. Possible origins of this effect are discussed. 
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INTRODUCTION 

Below the second roughening temperature 
(T < Tr 2 ~ 0.9if), the helium crystals having 
appeared in a metastable superfluid grow in two 
qualitatively different modes depending on initial su¬ 
persaturation Dpo 0, i- Until the deviation from 
equilibrium is rather low, the crystal face growth rate 
is low and determined by the well-known mechanisms, 
namely, growth on screw dislocations and Frank-Read 
sources [3|. However, the growth rate increases jumpwise 
above a certain critical supersaturation pc, which de¬ 
pends on temperature 0] , and a crystal grows completely 
in 200 — 400/is 0. This effect is characterized by a time 
lag between the nucleation of a crystal and a jumplike 
increase in its growth rate 0. The crystal grows at an 
ordinary rate in a rather long time, which also depends 
on the initial supersaturation and reaches 0.2s, and 
its growth then accelerates. At low temperatures, the 
growth rate increases by a factor of 100 — 1000. To date, 
the conditions under which the second growth scenario 
(burst-like growth) occurs have been determined 14 , the 
burstlike growth rates have been measured 0 0|, and 
the return to the normal state with the low growth rates 
typical of equilibrium crystals has been studied Q- 

However, the entire set of experimental data did 
not provide an unambiguous answer to the question of 
whether this phenomenon is the transition of a surface 
into a new rapidly growing state or is caused by a change 
in the volume properties of a crystal 0. Since the face 
growth kinetics changes, it is reasonable to assume that 
the effect is caused by a surface phase transition or the 
appearance of a qualitatively new intense growth source. 
As follows from filming the crystal growth at the stage of 
burstlike growth 0, the anisotropy of kinetic faceting 
is low, which means simultaneous acceleration of the 
growth kinetics of all faces (including various orienta¬ 
tions). This also means that the preparation for the be¬ 
ginning of burstlike growth should take the same time. 


Thus, the new growth mechanism should have a low sen¬ 
sitivity to the surface parameters such as the roughen¬ 
ing temperature, the surface energy, and the interplanar 
spacings. 

The difficulties in searching for such a mechanism led 
to assumption that the growth kinetic is accelerated due 
to a change in the internal state of the crystal. Only one 
experiment related to the volume parameters has been 
performed to date: researchers tried to detect a thermal 
effect during the formation of a rapidly growing state. At 
temperatures above 0.4 K, no specific features were re¬ 
vealed during a transition into the burstlike growth phase 
against the background of temperature oscillations ac¬ 
companying crystal growth. These results only give the 
upper estimate of the volume transition energy [8|. 

A change in the internal state of the crystal can af¬ 
fect its kinetic properties, such as thermal conductivity 
and internal friction. The intrinsic dissipation losses of 
an equilibrium crystal duriM vibrations at a frequency 
of 15-80 kHz are very high 0llOj. The maximum damp¬ 
ing decrement of helium is 0.5. The temperature depen¬ 
dences of the dynamic modulus and the decrement and 
the effect of ^He impurities, plastic deformation, and an¬ 
nealing on internal friction suggest that the dissipation 
mechanism has a dislocational nature 0- Since the in¬ 
ternal friction is determined by the sample volume, we 
can expect that the volume changes resulting in burst¬ 
like growth will affect the internal friction parameters of 
the crystal. Then, the return to the normal state with 
slow growth kinetics should be accompanied by the at¬ 
tenuation of the addition to the decrement in time. 

The purpose of this work is to measure the internal fric¬ 
tion at a frequency of about 75 kHz after crystal growth 
in the normal and anomalously fast (burstl ike) growth 
modes. Preliminary results are published in [ll| . 
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EXPERIMENTAL 


At the initial time, the liquid in a container was under 
the pressure exceeding the phase equilibrium pressure by 
DpQ. The value of Dpo was limited by spontaneous nucle- 
ation of a crystal on the wall, so that the maximum super¬ 
saturation was 5-12 mbar in various experiments. Before 
the onset of growth, the pumping of a quartz resonator 
was switched off and the amplitude of its vibrations be¬ 
gan to decrease. Then, a high voltage pulse was applied 
to a needle located so that the formed crystal touched 
a quartz oscillator. The additional electrostatic pressure 
stimulated the formation of a nucleation center for a solid 
phase 0, and the metastable liquid began to solidify. 
The crystals grown by this technique have a dislocation 
concentration of 10® — 10® cm~^ [3]. Crystal growth was 
filmed with a video camera. The switching-off signal from 
the generator was used to start the recording of quartz 
vibration amplitude and the readings of a capacitance 
pressure sensor at a step of 1 ps. The total recording 
time was 130 ms. The frequency and Q factor of the 
composite system were determined from 1-ms recording 
segments of vibration amplitude relaxation. The details 
of the optical technique were described in (o. [l2|. 


Figure 1 shows the quartz resonators used in exper¬ 
iments and the position of needles. The quartz fork 
had a resonance frequency of 31909.2 Hz and an intrin¬ 
sic damping decrement Sq « 10“'* in superfluid helium 
(Fig.la). Cylindrical quartz 3.4 mm in diameter and 
28 mm in length with a fundamental torsional vibration 
mode in helium was excited at a frequency of 74 551.9 Hz 
(Figs.lb,Ic). The Q factor of free vibrations was 3.2x 10*. 
The quartz vibration amplitude was chosen so that the 
maximum shear strain of a helium crystal did not ex¬ 
ceed 10“^. This resonator was used earlier to study the 
internal friction in crystalline helium i[i3. 


The helium single crystal having grown in the gap dur¬ 
ing torsional vibrations undergoes complex deformation, 
so that its torsional rigidity is determined by all elasticity 
tensor components. To process the experimental data, 
we use a simplified model: a helium crystal is consid¬ 
ered as an isotropic medium. The crystal is assumed to 
have the shape of a disk with its axis coinciding with the 
axis of torsional vibrations. In this approximation, the 
equation describing the vibrations of the quartz-helium 
crystal system has the form 



FIG. 1. Types of quartz resonators and the position of the 
tungsten tip initiating crystal growth. The arrows in (a) and 
(b) indicate crystal boundaries, (a) Crystal in the gap of a 
quartz fork, (b) crystal growth in a gap of 0.4 mm, and (c) 
assembly with central location of a tip and a gap of 0.18 mm 
(photograph and schematic section). 


volume Vc by the relationship 

Vc = Vo^kiDpo = Sh, (2) 

Ap 

where Vq is the inner container volume, p is the liquid 
helium density, Ap is the difference between the solid and 
liquid helium densities, fc; is the compressibility factor of 
liquid helium, and h is the gap size. Wavevectors ko and 
k of the torsional vibrations of quartz and helium are 
described by the conventional relations 


tan(A:oAo) / S' coi(kh) 

koLo kh ’ 

where Lq is the quartz crystal length and So and S are 
the quartz end face area and the helium crystal area, re¬ 
spectively. Helium crystal area S in contact with the end 
face after the end of growth was estimated from crystal 


fco = —, fc = —, c = 
Co c 


G' + iG" 
P' 


( 3 ) 


where cq and c are the torsional vibration velocities of 
quartz and the helium crystal, respectively; uj is the sys¬ 
tem vibration frequency; G = G' + iG" is the complex 
dynamic shear modulus of helium; and p' is the solid 
helium density. The imaginary part of the frequency 
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uj" = uj'512'k was determined from the experimental val¬ 
ues of the frequency and damping decrement of the sys¬ 
tem, and the complex frequency was then substituted 
into Eq. O- By solving the equation, we found the real 
and imaginary parts of the dynamic modulus of helium. 
The decrement of the medium was determined using the 
well-known relation 5He = ttG" jC. Using the simplifi¬ 
cations described above, this processing give results av¬ 
eraged over all elastic constants of the helium crystal. 

In earlier internal friction experiments [^, IT^ , we used 
a relation deduced from energy considerations to find 
5He- Under the conditions 5 27^ and Sq <SC S, the 

decrement of the medium is calculated by the formula 


^He 


S-5o 
1 - 


( 4 ) 


where ujq is the vibration frequency of free-standing 
quartz and a is the coefficient that takes into account the 
inertia of the helium crystal. This coefficient is a ~ 0.9 
for the given geometry and a gap of 0.4 mm. The decre¬ 
ments calculated by Eq.(|3]) and Eqs. © and m agree 
with each other to within 10%. 

The growth and the size oscillations of rapidly grow¬ 
ing crystals end in 0.5-3 ms depending on temperature 
Q. Then, the crystal volume increases due to the flow 
from an external system, and this increase is only 2% 
in 20 ms. Therefore, the contact area S is assumed to 
be constant in approximately 2 ms and can be calcu¬ 
lated by Eq.([2|). When calculating the modulus and the 
decrement of slowly growing crystals, we took into ac¬ 
count the change in their volume determined from the 
time dependence of the pressure decrement in the con¬ 
tainer. Contact area S{t) increased in time. Note that 
the estimation of S (t) by Eq. ([5]) gives a high error at the 
beginning of growth, when the crystal size is small or 
comparable with the gap size. 


EXPERIMENTAL RESULTS 
Experiments with the Quartz Fork 


In the experiments with the fork, spontaneous nucle- 
ation of the crystal on the container walls limited the 
starting supersaturation, which was at most 5 mbar. 
During crystal growth, the system frequency increased 
from 32 to 60 kHz. As follows from Eq.Q, the Q fac¬ 
tor of the system was mainly determined by the crystal 
with a decrement of 0.1-0.5. As a result, the vibrations 
of the fork-crystal system decay in 1-2 ms. In this time, 
the anomalous state has no time to form at such a low 
supersaturation Q. These experiments did not lead to 
the desired results. 


Experiments with Cylindrical Quartz 

In contrast to the experiments with the fork, the shift 
of the quartz-helium system frequency and the decay of 
quartz are low. Free vibrations of the system last the 
time that is sufficient for the anomalous state to form 
and the evolution of its internal friction to be studied 
during 130 ms. In the first series of experiments, a nee¬ 
dle was placed from one side and its tip was near the axis 
of quartz vibrations (Fig.lb). Experiments were carried 
out at temperatures of 0.49, 0.54, 0.59, 0.62, 0.65, 0.69, 
and 0.75 K to an initial supersaturation of 10 mbar. The 
scatter of the temperatures inside each group of mea¬ 
surements was 0.002 K. The gap between the end face of 
quartz and the metallic base was 0.4 mm. In the second 
series of experiments, a needle was located from below 
along the axis of torsional vibrations (Fig.lc). Experi¬ 
ments were performed at a temperature of 0.512 K, and 
the gap in this case was 0.18 mm. The position of the 
needle ensured the location of a crystal that corresponds 
better to the simplified model of calculating the crystal 
parameters. 

The averaged shear modulus calculated from the vol¬ 
ume compressibility and the Poisson ratio {v — 0.3) is 
1.22 X 10® dyn/cm^. The experimental values fall in the 
range (1 — 3.5) x 10® dyn/cm^. This difference in the 
results is related to the crystal anisotropy and the sim¬ 
plified approach to determining the shear modulus. In 
Eq.©, we did not take into account the presence of a 
rigid needle and the shift of the crystal from the axis of 
quartz vibrations. The last factor can be taken into ac¬ 
count by introducing a coefficient, which depends on the 
ratio of the crystal radius to the shift, into Eq.(IT]). This 
coefficient is unity for a zero shift. When the shift of the 
crystal is not taken into account, its torsional rigidity is 
underestimated and, hence, the shear modulus calculated 
by Eq. © is overestimated. The same effect is caused by 
a needle located from one side. This note belongs to the 
elastic modulus rather than to its time dependence, since 
the crystal size and shape change insignificantly in the 
measurement time. The decrement of a helium crystal 
calculated by Eq. © or Eq. o has a weak sensitivity to 
the shift. A numerical simulation showed that the shift 
of the crystal from the vibration axis even by its radius 
changes the decrement calculated by Eq.(|4|) by at most 
10 %. 


Slow Crystal Growth 

Figure 2 shows the changes in the pressure in the con¬ 
tainer and the frequency and decrement of the quartz- 
crystal system at T = 0.75 K and an initial supersatu¬ 
ration of 2.7 mbar. During crystal growth, the pressure 
in the container decreases in approximately 50 ms. The 
crystal radius calculated in approximation Q increases 
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FIG. 2. Time evolution of the quartz-helium crystal system 
parameters and the calculated values of (a) damping decre¬ 
ment and (b) shear modulus at T = 0.75 K in the normal 
growth region, (c) Change in the pressure in the container 
and the crystal radius calculated by Eq.([2]). (b) Recorded 

frequency of a quartz-helium crystal composite oscillator and 
the shear modulus calculated by Eq.((T]). (a) Change in the 
damping decrement of the system and calculated decrement 
of a helium crystal. 
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FIG. 3. Time evolution of (b) frequency and (a) damping 
decrement of the quartzhelium crystal system during crystal 
growth at various initial values of supersaturation for T = 
0.49 K after the end of burstlike growth. The numerals in 
(b) indicate the values of Dpo [mbar]. The sequence of the 
decrement curves from bottom to top in (a) corresponds to 
the sequence of the frequency curves in (b). 


and then, after a time of 20 ms, changes weakly. The 
frequency and decrement of the quartz-crystal system 
increase monotonically. The results of calculating the 
decrement and the shear modulus of a helium crystal are 
also depicted in Fig.2. It is seen that the shear modulus 
falls in the range (2.5 —3.5) x 10® dyn/cm^ during crystal 
growth. Allowing for the simplifications and the scatter 
of the experimental results, we assume that the shear 
modulus remains unchanged. The calculated decrement 
of a helium crystal changes weakly during crystal growth 
(Fig.2a). The absolute value of the decrement (0.2-0.3) 
agrees with the measurement results in [l^ . 

Thus, the slow crystal growth is characterized by a 
constant shear modulus and an unchanged, or weakly 
changed, decrement, which is close to the decrement of 
the crystals grown and annealed near the melting curve 
The growth from a metastable liquid weakly af¬ 
fected the internal friction of helium crystals. 


Burstlike Crystal Growth 

In the range of burstlike growth, the pressure detected 
by the pressure sensor decreases sharply after the nucle- 
ation of a crystal and exhibits a jumplike increase in the 
growth kinetics of all crystal faces. The decrement and 
frequency of the quartz-helium system were determined 
from recording free vibrations after completed crystal 
growth. Figure 3 shows the time dependences of these 
parameters. The following two regions are clearly distin¬ 
guished in these curves. After the end of crystal growth, 
the decrement of the system begins to decrease simulta¬ 
neously with an increase in the averaged dynamic shear 
modulus. In approximately 20 ms, relaxation is com¬ 
pleted and a further slow increase in the frequency and 
the shear modulus is observed. The second region is re¬ 
lated to crystal growth due to a liquid flow along a cap¬ 
illary from the external system. The slopes of the f{t) 
dependences in this region are proportional to initial su¬ 
persaturation Dpo, and the calculated shear modulus and 
the damping decrement of the crystal are constant within 
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FIG. 4. (a) Damping decrement and (b) shear modulus cal¬ 
culated from the curves in Fig. 3 and normalized by their 
steady values at t = 20 ms. The initial supersaturation is (□) 
8, (O) 5.7, (o) 4.5, and (A) 2.8 mbar. 

the limits of experimental error. The values of Sue in this 
region change from crystal to crystal and fall in the range 
0.05-0.25. For example, for the curves presented in Fig. 
3, these values are 0.074 (2.8), 0.1 (4.5), 0.13 (5.7), and 
0.11 (8); in parentheses, we give the values of Dp^ [mbar]. 
Such a scatter was observed earlier in annealed samples 
grown at a constant pressure and a constant temperature 
gradient, and it is likely to be related to different defect 
distributions in a sample [l^ . A nonuniform dislocation 
distribution in a crystal was detected earlier during slow 
growth, where the kinetics is controlled by defects (screw 
dislocations, Frank-Read sources Q.), from the deviation 
of a crystal shape from a regular hexagonal prism. The 
results obtained at the final stage of crystal growth agree 
with the well-known data. 

An excess decrement and a decrease in the dynamic 
modulus are only observed in the crystals grown in the 
burstlike growth mode. Figure 4 shows the evolution of 
the decrement and the shear modulus of a helium crystal 
normalized by their steady values at t > 20 ms. As 
is seen from Figs.3 and 4, the relaxation contribution 
increases with the initial supersaturation. The relaxation 
curves are well described by the exponential relationship 

= i + (5) 



FIG. 5. (a) Excess decrement amplitude and (b) process re¬ 
laxation time calculated by Eq.((5| vs. the initial supersatu¬ 
ration at a temperature of (□) 0.49, (O) 0.59, (o) 0.65, and 
(A) 0.69 K. 

The dependences of parameters A and r on Dp^ are pre¬ 
sented in Fig.5. The increase in the relaxation contri¬ 
bution is illustrated by an increase in amplitude A with 
the initial supersaturation. Parameter A vanishes at the 
supersaturation that corresponds to the boundary super¬ 
saturation separating the regions of normal and anoma¬ 
lous growth, A{Dpo —>• pc) ^0 0. The relaxation time 
falling in the range 2.5-4.5 ms depends on neither temper¬ 
ature nor the initial supersaturation within the scatter of 
the experimental results. 

Thus, the crystals grown in the burstlike growth mode 
exhibit an excess decrement and a decrease in the dy¬ 
namic shear modulus, which relax in time, in the tem¬ 
perature and supersaturation ranges under study. The 
effect is low or zero at the boundary between anomalous 
and normal growth and increases with the distance from 
it. This addition decreases at a time constant of 3 ms. 

DISCUSSION OF RESULTS 

The measured decay is caused by dissipation both in¬ 
side the crystal and at its surface. It is known from the 
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experimental data that high face kinetics after crystal 
growth lasts at least 3 ms and that the return to the 
normal slow kinetics occurs in 50-100 ms Q- The elas¬ 
tic stresses at the crystal boundary change the chemical 
potential of the solid phase and initiate interface motion. 
A similar effect was studied in [1^. As a result, melting 
takes place and crystal faces grow at a double frequency, 
since the addition to the chemical potential is quadratic 
in stress. In turn, the surface-assisted dissipation con¬ 
sists of two parts. The first part is caused by the surface 
growth kinetics, is determined by kinetic growth coeffi¬ 
cient AT, and is lower than 10”^. Oscillating growth is 
accompanied by liquid flows, which induce sound emis¬ 
sion at a double frequency. The second contribution is 
also low, less than 10“^. The sum of both contributions 
is much lower than the decrement of helium detected in 
experiments, which is at least 0.07. Thus, the additional 
decay is related to volume processes. 

The internal friction measured in experiments is deter¬ 
mined by the sum of all dissipation processes occurring in 
the crystal. The possible sources of the effect under study 
are as follows. The additional decay can be caused by the 
appearance of a new relaxation process. The time behav¬ 
ior of decaying agrees qualitatively with the change in the 
face growth kinetics after the end of growth. Therefore, 
the jumplike state of the crystal is assumed to manifest 
itself in both the appearance of an additional decay pro¬ 
cess and the acceleration of face growth. 

The second possibility consists in the effect of the 
changed state of the crystal on the oscillation of dislo¬ 
cation segments, which are considered to be responsible 
for the high damping decrement in crystalline helium. In 
the Granato-Liicke model of dislocational internal fric¬ 


tion [1^, the frequency and temperature dependences of 
internal friction in the case of damped oscillation of dis¬ 
location segments are expressed by the formula 




^'^dis 


1 + {iOTdis)'^' 


tms B{T)L'^ ( 6 ) 


where B is the damping constant, L is the average dislo¬ 
cation segment length, A is the dislocation concentration, 
n is the orientation factor, and ^ is the parameter that 
takes into account the length distribution of dislocation 
segments. For a 5 distribution, we have I; for a net¬ 
work of random intersections, ^ ~ 4.4. In terms of this 
model, a decrease in the decrement with a simultaneous 
increase in the dynamic modulus occurs when time Tdis 
increases, i.e., when the damping constant determined by 
the phonon gas viscosity increases. Whence, it follows 
that the dislocation motion in the crystal grown in the 
burstlike growth mode is retarded by phonons to a lesser 
extent as compared to that in an equilibrium crystal. 

The last conclusion holds true if the geometric con¬ 
figuration of a dislocation network remains unchanged 
after the end of growth. However, as is known from 
earlier experiments on internal friction [l^ . even slow 


growth of helium crystals can lead to time variations of 
the decrement and the dynamic elastic modulus of a crys¬ 
tal. It is necessary to perform long-term annealing of a 
helium crystal near the melting curve and to provide a 
slow temperature change during measurements to obtain 
reproducible results. Taking this fact into account, we 
assume that burstlike crystal growth creates a nonequi¬ 
librium network of dislocations and it begins to relax af¬ 
ter the end of growth, which is reflected on the values of 
the decrement and the dynamic modulus. For a network 
of dislocations with self-intersections, we have AL^ « I; 
therefore, it weakly changes upon annealing and the time 
evolution is related to an increase in the average dislo¬ 
cation segment length during rapid annealing, according 
to Eq. ([6|) . This assumption is in conflict with the results 
of burstlike growth experiments on the crystals that have 
no Frank-Read sources and screw dislocations [3,0. In 
these experiments, a crystal face had no growth sources 
after the end of burstlike growth; that is, the burstlike 
growth mode did not generate defects. Even multiple 
repetition of the process did not degrade the surface. 

In conclusion, note that the excess decaying and the 
modulus defect, which relax in time, are only observed 
in the crystals grown in the burstlike growth mode. It 
is impossible to find the source of the anomaly using the 
data obtained. The effect can be explained by both a 
change in the internal state of the crystal, which mani¬ 
fests itself in either the appearance of a new relaxation 
channel or an action on dislocation drag, and annealing 
of a nonequilibrium network of dislocations, which ap¬ 
pears during burstlike growth. The question of whether 
burstlike growth can generate a nonequilibrium network 
of dislocations remains open. 
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